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Abstract 



We present a model of the lepton masses and flavor mixing based on the 
discrete group 5 4 ® Z 2 ■ In this model, all the charged leptons and neutrinos 
are assigned to the 3 a representation of in the Yamanouchi bases. The 
charged lepton and neutrino masses are mainly determined by the vacuum 
expectation value structures of the Higgs fields. A nearly tri-bimaximal lepton 
flavor mixing pattern, which is in agreement with the current experimental 
results, can be accommodated in our model. The neutrino mass spectrum 
takes the nearly degenerate pattern, and thus can be well tested in the future 
precise experiments. 
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I. INTRODUCTION 



To understand the origin of fermion masses and flavor mixing is crucial and essential in 
modern particle physics. From the analyses [1] of recent neutrino oscillation experiments 
[2-6], we have confirmed the large solar mixing angle 9 12 ~ 34° [7], maximal atmospheric 
mixing angle 9 23 — 45° and a very tiny 6 13 , l3 < 10°. A very small solar mass squared 
difference Am^ = m 2 — m 2 ~ 8.0 x 1CT 5 eV 2 and an atmospheric mass squared difference 
Am| 2 = m\ — wij ~ ±2.5 x 10~ 3 eV 2 are given by the experimental data [1]. Here the plus 
and minus signs in front of Am| 2 correspond to the normal mass hierarchy (m 1 < m 2 < m 3 ) 
and inverted mass hierarchy (m 3 < m 1 < m 2 ) case. Since the standard model predicts 
massless neutrinos, it is clear that we should extend the standard model to accommodate 
non-vanishing neutrino masses. Among all possible mechanisms, the seesaw mechanism [8] 
is a very interesting and elegant one to explain the light neutrino masses. In the framework 
of the seesaw mechanism, the light left-handed neutrino masses can well be understood by 
introducing heavy right-handed Majorana neutrinos, and the light neutrino mass matrix is 
given by 

M v = M D M^M T D , (1) 

where M D and M R are the Dirac and Majorana mass matrices respectively. The typical 
mass scales of M R are 10 14 ~ 10 16 GeV. 

Since the Yukawa structures are not well constrained up to now, to identify the structure 
of the neutrino mass matrix is one of the main objects in neutrino physics. An interesting 
and natural way to study the Yukawa coupling and find the underlying physics is the flavor 
symmetry. Among all possible flavor symmetries, discrete non-Abelian groups such as S 3 
[9], A A [10] and so on have attracted a lot of attention. Note that most of these models rely 
on the basis of their group representation, and contain more theoretical parameters than 
the observables. 

The permutation group S 4 , which is formed by the 4! permutations, totally contains 
24 group elements which belong to five conjugate classes. Therefore it has five irreducible 
representations (reps). Among these irreducible reps, there are two one-dimensional (l s and 
1 A ), one two-dimensional (2) and two three-dimensional (3^ and 3g) reps. Here the 'S' and 
'A' mean symmetric and antisymmetric reps respectively. The character table of S 4 is given 
in Table 1. In respect that the particles in our model are all assigned to the l s , 2 and 3 a , 
we list the relevant representation matrices in Appendix A. 

The S A models for fermion masses have been discussed by several authors in Refs. [11,12]. 
In this note, we adopt the S 4 <8> Z 2 flavor symmetry and assign all the charged leptons, 
neutrinos and Higgs into the irreducible Yamanouchi bases [13]. Such a compact scheme 
contains only a few parameters, and therefore it can be examined quite well in the future 
experiments. 

In the following, we will present the main contents of our model in section II. Some 
detailed analytical and numerical analyses will be given in section III. In section IV, the S4 
invariant Higgs potential are discussed. Finally, a brief summary is given in section V. 
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II. PARTICLE ASSIGNMENT 



In our model, all the leptons are assigned to the 3^ rep of S 4 , but their Z 2 charges are 
different. Under the group 

G = SU(3) C <g> SU(2) L <g> U(1) Y ®S a ®Z 2) (2) 

the lepton contents in our model are placed as 

€ L ~(l,2,-l)(3 a ,+), 
e R ~(l,l,-2)(3 a ,+), 

^-(l^O)^,-) , (3) 

where the plus or minus sign denotes the reflection properties under Z 2 , i.e. ^ — > The 
Higgs scalars in our model are placed as follows 

tf u -(1,2,-1)03,+) , 
ff d ~(l,2,-l)(lg,-) , 

x-a^-ix^,-) , 

0~(1,2,-1)(2,+) . (4) 

Here the Z2 symmetry guarantees that the lepton doublet couples to (e#, fa) and (z/r, %j) 
respectively. It will be shown in section IV that the masses of Higgs scalars are almost 
unconstrained in this model, thus we can choose their masses at some high energy scales in 
order to avoid the tree level flavor charging neutral currents. Note that we do not introduce 
any SU(2) L singlet or triplet Higgs [14]. Hence the Heavy right-handed neutrino masses are 
exactly degenerate. 

By using the group algebra given in Appendix A, we write down the S 4 <E> Z 2 invariant 
Yukawa couplings: 

-C Y = a e £ iL e iR H u + PJ ijk £ iL e jR (p k + aJ iL v iR H d + P u g ijk £ iL v jR x k + \ M ^R v iR + h.c. (5) 

where H u = iT 2 H* and <fi = ir 2 <f)*. The last term in Eq. (5) is the bare Majorana mass term 
with M being the typical mass scale of the heavy right-handed neutrinos. The coefficients 
a e v and (3 are in general complex parameters. However, since their phases are all global 
phases, there will be no Dirac CP violating phases in the MNS matrix and only the Majorana 
phases can be accommodated up to now. 1 The structures of the traceless matrices / and 
g can be obtained from the CG coefficients of S 4 . Following the CG coefficient tables given 
in Ref. [13], we arrive at 





/«i = V2 1 , f ij2 = -1 , (6) 



Here we do not consider the spontaneous CP violation [15] 
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V-l -y/2 / 

where all the coefficient matrices are symmetric and traceless. Hence both the charged 
lepton and neutrino mass matrices are symmetric. 

As we have mentioned above, the right-handed neutrino masses are exactly degenerate 
by now. Such a result is not favored in some leptogenesis models [16] in which the mass 
splits of the heavy right-handed neutrinos are required. This drawback can easily be solved 
by introducing an additional Higgs scalar rj ~ (1, 1, 0)(2, +) or £ ~ (1, 1, 0)(3 a , +) which 
couples to u R in the form of 

~ C rt-u R = PllfijkVmVjRVk , ( 8 ) 

or 



- C Z-v R = ^R9ijk v iR v jRik ■ (9) 

Their vacuum expectation values(VEVs) will bring the mass splits of right-handed neutrinos 
and the phases of f3 R will lead to the Dirac CP violating phases in the MNS matrix. 

Another possible way to gain the mass differences of the right-handed neutrinos is to 
consider the renormalization group equation running effects [17] which can usually be used 
to generate a small mass split in the resonant leptogenesis models [18]. A Dirac CP violating 
phase can also be gained at the same time from the renormalization group equation running 
[19]. 

In our analyses, we focus our attention on the low energy phenomena of our model. So 
we assume that all the Yukawa couplings are real and the right-handed neutrino masses are 
exactly degenerate. 



III. LEPTON MASS MATRICES AND MIXING MATRIX 



Assuming the flavor symmetry S 4 is broken by the VEVs of the Higgs scalars in the form 
of (H u ) = v u , (H d ) = v d , = u { and (xi) = w { , we obtain the mass matrix of the charged 
leptons as 



( a 
V / 2x 1 a + x 1 



.1' ■-) 



x 2 
a — x. 



(10) 



where a and x { are defined by a = a e v u and x { = f3 e u i respectively. 
The Dirac mass matrix of the neutrinos is given by 
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M D = 



/2 yi + b -y 2 -y 3 \ 

-2/ 2 V2y 2 ~Vi + b ~V2y 3 

V -y 3 -V2y 3 -Vi - V2y 2 + &/ 



(11) 



where b = a v v d and y i = P^w^ Taking account of the seesaw relation in Eq. (1), we obtain 
the light neutrino mass matrix: 



M v = m Q 



({2 + e) 2 + r 2 + s 2 V2s 2 - (1 + 2e + y/2r)r -(1 + 2e - 2y/2r)s \ 

(1 -e- v / 2r) 2 + r 2 + 2s 2 (2y/2 - 2^2e + r)s 

(1 -e + v / 2r) 2 + 3s 2 ; 



(12) 



where m = yf/M, e = b/y 1 , r = y 2 /y 1 and s = y 3 /y v Note that all the lepton mass 
matrices are symmetric. 

The mass matrices M t and M v can respectively be diagonalized by two unitary matrices 
V l and V v : V l M^* = Diag{m e , m^, m T }, VjM u V* = Diag{m 1 ,m 2 ,m 3 }. The lepton flavor 
mixing matrix Vm NS , which links the neutrino mass eigenstates u 2 , v> 3 ) with their flavor 
eigenstates (y e , i/ T ), is then given by V MNS = V}V V . 

A. charged lepton sector 

By diagonalizing the charged lepton mass matrix in Eq. (10), we get the charged lepton 
masses: 



m e = a + 2x l , 

m = a — x 1 + v^3x 2 , 



\/3. 



x. 



2 5 



(13) 



and the mixing matrix Vf 



2 

v 



1 


1 \ 






1 


1 




"7e 


i 


1 


71 


"72/ 



(14) 



A typical interesting and instructive feature is that V l is a constant matrix with respective 
to arbitrary a and x { . After being transposed, this matrix takes the very similar form to the 
tri-bimaximal mixing matrix [20] 



^tri-bi 



( - 
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y/2/ 



(15) 
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only up to an exchange between the first and second column. 

By using Eq. (13), we can also express the parameters a, x 1 and x 2 in terms of the 
charged lepton masses as 



a = ^ ( m e + "V + m r) , 

% = ^= («V - «v) • (16) 

Considering the mass hierarchies of the charged leptons: m T ^> ^> m e , it is easy to 
conclude that \a\ > \x 2 \ > \x ± \. 



B. neutrino sector 



As an example, we assume s = and r < 1 2 . Then, M D takes the simple form: 

'2 + e -r \ 

Afr, = 2/i | -r -l + V2r + e 

-1 - ^2r + e / 



(17) 



and the light neutrino mass matrix can be written as 



[2 + ef + r 2 -r- 2re - y/2r 2 
•-2re- v / 2r 2 (1 - e - V2r) 2 + r 
















(1 - e + V2r 



(18) 



There are two texture zeros in the third row of M v [21]. Therefore only the mixing between 
the first two generations exists in the neutrino sector. After straightforward calculations, 
we obtain the light neutrino masses: 



m 1 ~ m 



m 2 ~ m 



m 3 ~ m 



(l-e) 2 -2v / 2(l-e)r + -(4-e) 



(2 + 



■(2 + e)r 2 



(1 - e) 2 + 2V2(1 - e)r + 2r 2 



(19) 



where the terms in proportion to r 3 have been neglected. The mixing matrix V u can be 
approximated to 



2 In the case of s ^ 0, there will be more freedom for us to fit the data in principle. However, it will 
involve the problem of cubic root. Since the main propose of this paper is to show the feasibility 
of this 5 4 model, we simply take s = as an example. In our numerical calculations at the end of 
section III, we do not make any approximation on r or s, and we will see that the exact numerical 
result satisfies the condition s,r«l very well. 
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r v — 
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-r H r 
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1 2 \ 

-1 + — r 2 
18 

1 ^2 2 

-r H r 

3 9 

1/ 



+ 0(r 3 ) . 



(20) 



It is interesting that only the parameter r appears in the approximate expression of V v . 
The MNS mixing matrix is given by 

1 0' 

1 
1 

\ 



^MNS = V}V V 










1/ 
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where 



sin i 



V2 



cos 9 = 1 r 
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(21) 

(22) 
(23) 



The last matrix in Eq. (21) contributes a trivial phase ±ir in the MNS matrix, thus it can 
be safely neglected without changing any phisics. The explicit expression of V^ MNS is 



2 cos 9 sin 9 



^MNS — 



■ cos 9 
cos 9 



sin 9 
sin 9 



cos 9 2 sin 9 
+ 



cos 9 
cos 9 



V6 
sin 9 

sin^ 



V 



Vq Vz Vz Vq 

Comparing Eq. (24) with the standard parametrization 
( 

^MNS = 

V S 12 C 23 ~~ C 12 C 23 S 13 e 

we obtain three mixing angles 



X 
1 

1 

"71/ 



(24) 



c 12 c 13 
— S 12 C 23 ~ C 12 S 23 S 13 e 



S 12 C 13 

iS 

C 12 C 23 ~~ S 12 S 23' S 13 e 
~ C 12 S 23 — S 12 C 23 S 13 e 



s 13 e- iS \ 


/ e i P 





0\ 


S 23 C 13 | 


■ 







C 23 C 13 / 


V o 





1/ 



(25) 



sin# 



cos# 2sin# 



12 



V3 
1 



S ' m9 23 = ^ ' 

sin 13 = , 



(26) 
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together with 5 = p = o = 0. 

From the experimental data of 8 l2 and Eq. (26), we can estimate that 8 ~ —1.4°. These 
results depend on ours assumptions. It should be emphasized that 8 13 will not vanish if we 
loose our assumptions by setting s ^ 0. In this case, small corrections on 8 23 will also be 
given corresponding to the nonzero s, and we will illustrate these corrections in the next 
subsection. 



C. numerical analyses 

In our numerical analyses, we first take the charged lepton pole masses Mf olQ given by 
PDG [22], and then evaluate the charged lepton running masses m^/i) from Mf olc to the 
M z scale by using [23] 



mi (M z ) = Mr le 

Then we obtain 



vr \2 mi {M z ) / 



(27) 



m e (M z ) = 0.4867 MeV , 
m^M z ) = 102.7 MeV , 

m T (M z ) = 1746.6 MeV . (28) 
We can directly calculate the parameters (a, £1,2:2) by using Eqs. (16), 

a ~ 616.6 MeV , x 1 ~ -308.0 MeV , x 2 ~ -474.5 MeV . (29) 



For the neutrino sector, we take the light neutrino mass squared differences Am 



2 

21 



8.0 x 10~ 5 eV 2 and Am 2 2 = ±2.5 x 10~ 3 eV 2 , three mixing angles 12 = 33.9°, 6 23 = 45° 
and 8 13 = 0° [1] as the input values. We firstly consider the case s = 0. By combining Eqs. 
(19), (22) and (23), after some numerical calculations, we find that only the normal mass 
hierarchy case is allowed for our assumptions. The corresponding parameters are 



m 



~ 2.99 x 10~ 2 eV , r ~ 6.92 x 10~ 2 , e ~ -0.545 . (30) 



It can be seen clearly that our assumption r <^ 1 is quite reasonable. The lightest neutrino 
mass is m 1 ~ 0.063 eV. Hence the neutrinos have a nearly degenerate mass spectrum: 
m 1 ~ m 2 ~ m 3 . As for the case s ^ 0, the mixing angles #13 and 823 will in general get 
corrections and deviate from their tribi-maximal values. In Fig. 1, we show the corrections 
on #13 and 823 from nonzero s. We can find that 813 is more sensitive to the non-vanishing 
s, and in order to get a nonzero 813 within 1 and 1o confidence level s should not be larger 
than 0.02 and 0.035. 

Note that the predictions of our model rely on the VEV structures of the Higgs fields. 
Here we have assumed that such conditions on the VEVs of the Higgs scalars can be well 
satisfied. In section IV, a detailed analysis will be given to show that this kind of VEV 
configuration can be obtained by a suitable choice of parameters of the Higgs potential. As 
mentioned above, in our analyses we do not consider the CP violation of our model. To 
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generalize our model to the CP violating case, we can add the singlet Higgs scalar given 
by Eq. (8) or (9) and insert the phase factors into the parameters a and (3 in Eq. (5). In 
this case, the mass splits of the right-handed neutrinos and the CP violation in the neutrino 
mixing can be acquired simultaneously. As a consequence, there will be more freedom to fit 
the experimental data and one may test the leptogenesis mechanism within this framework. 
Such an extension may be more interesting and the corresponding detailed analyses will be 
elaborated elsewhere. 



IV. HIGGS POTENTIAL 



The most general S4 <g> Z 2 invariant Higgs potential in our model is given by 
V= -ml(HlH u ) - m d (H d H d ) - mj £ - m 2 ^ X )xj 



i=i 
2 



j'=i 



+ \ u {HlH u ) 2 + \ d {HlH d f + 4^ + \ x (^2 x)xj 

+ Ti [((/>[(/>! - 4<f>2) 2 + (0l02 + 4<Pl) 2 } + T 2 (4<f>2 - 020l) 2 

+ A*i [( 2 XiXi - X2X2 - xlxs) 2 + (X1X2 + X2X1 - V2X2X2 + V2X3X3) 2 

+ (xlxs + X3X1 + V2xlxs + V2xtx2) 2 } 

+ V2 [(X1X2 - X2X1) 2 + (XiXs - x\x\) 2 + (xlxs - X3X2) 2 

+ V3 [(X2X2 - X3X3 + V2X1X2 + V^xixi) 2 + (X2X3 + X3X2 - V2x\x3 - V2xlxi) 2 } 
+ XMH^HlH,) + \ 2 {HlH d ){H\H u ) + [UHlH d f + h.c] 



+ \t(HtH u ) (X>&J + {4 [{Hl^f + [Hl<p 2 f] + h.c.} + \l (|/Jt 01 | + \Hl<p 2 

+ {\t [(Hl^M^ - <&<j> 2 ) - (Hl<h)(<l>[<h + 4<f>i)} + h.c.} 

+ ^{HlH u ) (Ex!*) + { A 2 [(^Xi) 2 + (#k 2 ) 2 + (Hlx,) 2 ] + h.c.} 

+ A^(|iJt Xl | 2 + | jf 7t X2 | 2 + | iJ t X3 |^ 

+ <rf (|>&) + H [(^]0i) 2 + (^]0 2 ) 2 ] +h.c.} + af(|ijt 01 | 2 + #t 
+ (l>b) +{^2 X [(^Xi) 2 + (^X2) 2 + (^X3) 2 ] +h.c.} 

+ a* (|^xi| 2 + |^k 2 | 2 + |^3x3 

+ {^4 [(^]xi)(2XiXi - X2X2 - X3X3) - (#]x2)(XiX2 + X2X1 - ^2X2 + V2X3X3) 
- {H\ X 3){x\x3 + X3X1 + V2 X lx3 + V2xtx2)} + h.c.} 



d V2 
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+ «i E 



vi=l 




Vi =1 / 

+ (0102 + 020l)(X2X3 + X3X2 - V2XlX3 ~ ^xtxi) 

+ {^3 [2(01x2 + 0ks) 2 + (01X2 - 4X3 + V2<p\ X i) 2 + (4X2 + 0^3 - v^Xi) 2 ] + h.c.} 

+ {^4 [(01X2 - 4X3 - V24xiY + (02X2 + 01X3 + ^Ixif + 2(4x2 " 01x3)'] + h.C.} 



402)(X2X2 - XsXs + v / 2x t iX2 + v^xki) 
h.c.} 



+ «5 
+ «6 



^lX2 + 02X3 2 + 0lX2 - 02X3 + ^01X1 2 + 02X2 + 0lX3 - V^Xl 



+ 2 



lX2 - 02X3 - V24lXl + 02X2 + 0lX3 + ^02X1 

+ [(#t0i)(x2X2 - xks + V2 X \x2 + v^xki) 

" (Hl<h)(x2X3 + X3X2 - v^xlxs - v^XsXi)] + h.c.} 

+ {5 2 [v / 2(^Xl)(0 t lX2 + 02X3) + (^X2)(0IX2 - + V2<j>\ X l) 

- (HtX3)(4lX3 + 02X2 - v^Xi)] + h.c.} 

+ {5 3 [V2(Ht Xl )(xi(f>i + X3> 2 ) + (^X 2 )(xl0i - X 3 02 + V2 X \<Pi) 

- (Htx3)(xl<i>i + xl<h-V2x[<h)]+h.c.} . 



02X2 - 0lX3 



(31) 



The Higgs potential contains totally 39 parameters. Among these parameters 12 of them 
are in general complex and the rest are real. Here we take all the parameters and VEVs to 
be real and the minimum of Higgs potential can be written as 



t / 2 2 2 2 2 / 2 , 2\ 2 / 2 , 2 , 2\ 

Knin = ~m u v u - m d v d - + u 2 ) - m x (w 1 + w 2 + w 3 ) 

+ Kvi + W d + (A + ri)K + u 4 2 ) + (A x + + (X x + 3^ + fi 3 )(w 4 2 + wf) 

+ . . . . (32) 

The full analytical formulae of V m i n have been listed in Table 2. Note that the couplings 
r 2 and fi 2 do not appear in the Higgs minimal. The number of parameters in the Higgs 
potential is more than 30, thus we are confident that it is adequate to arrive at the suitable 
minimum of the Higgs potential, and the VEV structure in section III can be satisfied. 



V. SUMMARY 

We have presented a lepton mass model based on the discrete group S 4 <g) Z 2 . The 
right-handed charged leptons and SU(2) L doublets are all assigned to the 3^ reps of S 4 
with plus Z 2 charges, and the heavy right-handed neutrinos are embedded in the 3 a with 
minus Z 2 charges. Four sets of Higgs doublets are introduced and the lepton masses are 
mainly determined by the VEV structures of the Higgs fields. After diagonalizing the lepton 
mass matrices, we obtained a nearly tri-bimaximal mixing matrix together with a nearly 
degenerate light neutrino mass spectrum. By some careful analytical and numerical analyses, 
we show that our model can well fit the current experimental values of three charged lepton 
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masses, two neutrino mass squared differences and three mixing angles, although it contains 
only a few free parameters. 

In some recent papers [12], the flavor S A symmetry has been treated as a subgroup of 
the continuous flavor 5*0(3) or SU(3) symmetry together with the gauge symmetry in some 
grand unified theories, i.e. the SO (10) models, and the quark mixing can be included. 
A simple way to contain the quark mixing in our model is to assume that all the quarks 
belong to the identity reps, and then the quark mixing can be obtained via the standard 
way. A supersymmetric extension of our model should also be interesting and may be given 
elsewhere. 
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APPENDIX A: REPRESENTATION MATRICES OF S A 



The 2 reps for the class C 2 in the Yamanouchi bases are 





(12) ~ (34) 

(13) ~ (24) 



-1 



(23) ~ (14) 



_ i 

2 
2 



2 



_1 

2 

2 



4' 
2 ■ 



(Al) 



where the definition of C 2 can be found in Table 1. All the other group elements can be 
obtained by using the relation 



(hh ' " h ' ' ' in-iK) — (hh ' " h) (h ' ' ' K-iK) 
The matrices of the 3^ reps are given by 



(12) 



(34) ~ & 



(14) 




(23) 



(13) 



i 

1 

2 




and the representation matrices in the 3^ are 



(12) 



(34) 



(14) 




1 


V3 


2 


2 


^3 


1 


2 


2 








/ " 


1 

2 




V3 




2 


V 






(24) 



















2 


1 
2 


■) 







1 


1 


2 


2V3 


1 


5 


2^/3 


6 




^2 




3 



/ —I 1 



1 

2^ 
1 
2 



'2\ 
3 > 

yj 

3 

1 

3 



(A2) 



(A3) 



(A4) 



The other group elements can be obtained by using Eq. (A2). 



12 



REFERENCES 



[1] A. Strumia and F. Vissani, hep-ph/0606054. 

[2] SNO Collaboration, Q.R. Ahmad et al, Phys. Rev. Lett. 89, 011301 (2002). 

[3] For a review, see: C.K. Jung et al, Ann. Rev. Nucl. Part. Sci. 51, 451 (2001). 

[4] KamLAND Collaboration, K. Eguchi et al, Phys. Rev. Lett. 90, 021802 (2003). 

[5] K2K Collaboration, M.H. Ahn et al, Phys. Rev. Lett. 90, 041801 (2003). 

[6] CHOOZ Collaboration, M. Apollonio et al, Phys. Lett. B 420, 397 (1998); Palo Verde 

Collaboration, F. Boehm et al, Phys. Rev. Lett. 84, 3764 (2000). 
[7] L. Wolfenstein, Phys. Rev. D 17, 2369 (1978); S.P Mikheyev and A.Yu. Smirnov, Yad. 

Fiz. Sov. J. Nucl. Phys. 42, 1441 (1985). 
[8] P. Minkowski, Phys. Lett. B 67, 421 (1977); T. Yanagida, in Proceedings of the Workshop 
on Unified Theory and the Baryon Number of the Universe, edited by O. Sawada and 
A. Sugamoto (KEK, Tsukuba, 1979), p. 95; M. Gell-Mann, P. Ramond, and R. Slansky, 
in Super gravity, edited by F. van Nieuwenhuizen and D. Freedman (North Holland, 
Amsterdam, 1979), p. 315; S.L. Glashow, in Quarks and Leptons, edited by M. Levy et 
al. (Plenum, New York, 1980), p. 707; R.N. Mohapatra and G. Senjanovig, Phys. Rev. 
Lett. 44, 912 (1980). 
[9] P. F. Harrison and W. G. Scott, Phys. Lett. B 557, 76 (2003). 
[10] E. Ma and G. Rajasekaran, Phys. Rev. D 64, 113012 (2001); K.S. Babu, E. Ma and 

J.W.F. Valle, Phys. Lett. B 552, 207 (2003). 
[11] S. Pakvasa and H. Sugawara, Phys. Lett. B 82, 105 (1979); E. Derman and H.S. Tsao, 
Phys. Rev. D 20, 1207 (1979); S. Pakvasa, H. Sugawara and Y. Yamanaka, Phys. Rev. 
D 25, 1985 (1982); R.N. Mohapatra, M.K. Parida and G. Rajasekaran, Phys. Rev. D 
69, 053007 (2004); E. Ma, Phys. Lett. B 632, 352 (2006). 
[12] C. Hagedorn, M. Lindner and R.N. Mohapatra, JHEP 0606, 042 (2006); Y. Cai and 

H.B. Yu, hep-ph/0608022. 
[13] J.Q. Chen, J.L. Ping and F. Wang, Group Representation Theory for Physicists, 2nd 

edition, World Scientific (2002). 
[14] W. Chao and H. Zhang, Phys. Rev. D 75, 033003 (2007). 
[15] T.D. Lee, Phys. Rev. D 8, 1226 (1973). 

[16] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986); P. Langacker, R.D. Peccei 
and T. Yanagida, Mod. Phys. Lett. A 1, 541 (1986); M.A. Luty, Phys. Rev. D 45, 455 
(1992). 

[17] PH. Chankowski and Z. Pluciennik, Phys. Lett. B 316, 312 (1993); K.S. Babu, C.N. 
Leung and J. Pantaleone, Phys. Lett. B 319, 191 (1993); S. Antusch, M. Drees, J. 
Kersten, M. Lindner and M. Ratz, Phys. Lett. B 519, 238 (2001); Phys. Lett. B 525, 
130 (2002). 

[18] M. Flanz, E.A. Paschos, U. Sarkar and J. Weiss, Phys. Lett. B 389, 693 (1996); L. 

Covi and E. Roulet, Phys. Lett. B 399, 113 (1997); A. Pilaftsis, Nucl. Phys. B 504, 61 

(1997); Phys. Rev. D 56, 5431 (1997). 
[19] S. Luo, J.W. Mei and Z.Z. Xing, Phys. Rev. D 72, 053014 (2005); Z.Z. Xing and H. 

Zhang, hep-ph/0601106. 
[20] PF. Harrison, D.H. Perkins, and W.G. Scott, Phys. Lett. B 530, 167 (2002); Z.Z. 

Xing, Phys. Lett. B 533, 85 (2002); PF. Harrison and W.G. Scott, Phys. Lett. B 535, 



13 



163 (2002); X.G. He and A. Zee, Phys. Lett. B 560, 87 (2003); C.I. Low and R.R. 
Volkas, Phys. Rev. D 68, 033007 (2003); E. Ma, Phys. Lett. B 583, 157 (2004); hep- 
ph/0409075; G. Altarelli and F. Feruglio, Nucl. Phys. B 720, 64 (2005); F. Plentinger 
and W. Rodejohann, Phys. Lett. B 625, 264 (2005); K.S. Babu and X.G. He, hep- 
ph/0507217; A. Zee, Phys. Lett. B 630, 58 (2005); S.K. Kang, Z.Z. Xing and S. Zhou, 
Phys. Rev. D 73, 013001 (2006); X.G. He, Y.Y. Keum and R.R. Volkas, JHEP 0604, 
039 (2006); N. Haba, A. Watanabe and K. Yoshioka, Phys. Rev. Lett. 97, 041601 (2006); 
Z.Z. Xing, H. Zhang and S. Zhou, Phys. Lett. B 641, 189 (2006); G. Altarelli, F. Feruglio 
and Y. Lin, Nucl. Phys. B 775, 31 (2007). 

[21] Z.Z. Xing and H. Zhang, Phys. Lett. B 569, 30 (2003). 

[22] W.M. Yao et al, Particle Data Group, J. Phys. G 33, 1 (2006). 

[23] H. Arason et al, Phys. Rev. D 46, 3945 (1992). 



14 



TABLES 



TABLE I. The character table of S* 4 . Here N stands for the number of elements in class C i 
and C denotes the cycle structure of each class. 



Class 


N 


C 


Is 


1a 


2 


3 a 




c, 


1 


l 4 


1 


1 


2 


3 


3 


c 2 


6 


21 2 


1 


-1 





1 


-1 


c 3 


8 


31 


1 


1 


-1 








c 4 


6 


4 


1 


-1 





-1 


1 


c 5 


3 


2 2 


1 


1 


2 


-1 


-1 



TABLE II. The Higgs potential at its minimum. The parameters Ki(i = 1 — 4) are de- 
fined as K\ = m + 4^3 + 4«4 + 2k^ + 2kq, K 2 = k\ — 2k 2 + + 6K4 + K5 + 3^6, 
K 3 = ki + 2k 2 + 6K3 + 2/t4 + 3«5 + K6 and K4 = 2k 2 + 2k 3 — 2k 4 + k§ — kq. 



Terms 


Coefficients 


Terms 


Coefficients 


v 2 

u 


-ml 




-m 2 


u 2 + u 2 


-ml 


tof + W 2 + U> 3 


-m\ 




A u 


v d 


Xd 


(u 2 + u 2 ) 2 


+Tl 


(w\ + W2 + w 2 ) 2 


Ax 


w 2 


4/ii 


4 1 4 
^2+^3 


3/il + ^.3 


W 2 W 2 + W 2 W 2 


8/i 3 


^2^3 


6/ii + 2fi 3 


W1W2 


4V2(Ai 3 - Mi) 


W\W 2 W 2 


12v^(a*i - Ma) 




Af + 2\t + Xf 




af + 24 + a* 


v 2 u {wl + wl + wl) 


Xf + 2A* + A* 


V 2 d {wi +W 2 2 + wl) 


fjf + 2cr£ + erf 




Ai + A 2 + 2A 3 


v u u\ 


2A? 


VuUlU 2 . 


-6Af 


v d wf 


4<T? 


v d wl 


2^2^ 


v d wi{wl + wl) 


-6a 4 x 


v d w 2 wl 


-6V2a* 


(uj + u 2 2 )wl 


K x 


u 2 w 2 + u 2 w 2 


K 2 


v\w\ + ^2^3 


K 3 


u\w\w 2 


2V2K 4 


U\U 2 W\W 3 


-4V2K 4 


U\U 2 W 2 W 3 


AK A 


U2W\W 2 


-2V2K4 


V u UlWlW 2 


4V2(<5i + <5 2 + <*s) 


v u u\w\ 


2(5 l +5 2 + S 3 ) 


V u Ulwl 


-2(5! + 6 2 + 5 3 ) 


V u U 2 WlW 3 


4^2(^1 + <5 2 + 5 3 ) 


V U U 2 W 2 W 3 


-4(<5i + 5 2 + <5 3 ) 
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FIGURES 




FIG. 1. Relations between nonzero s and mixing angles (#13, #23)- The dashed lines cor- 
respond to the 1, 2 and 3a confidence level from oscillation experiments. We also take 
Am§ 2 = 2.5 x KT 3 eV 2 , Am| x = 8.0 x lCT 5 eV 2 and 6 12 = 33.9° in our calculations. 
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